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Abstract 

We study the conformal bootstrap for a 4-point function of fermions {'ip'f'fJ') ™ 3D. 
We first introduce an embedding formalism for 3D spinors and compute the conformal 
blocks appearing in fermion 4-point functions. Using these results, we find general 
bounds on the dimensions of operators appearing in the ■0 x ■0 OPE, and also on 
the central charge Ct- We observe features in our bounds that coincide with scaling 
dimensions in the Gross-Neveu models at large N. We also speculate that other 
features could coincide with a fermionic GET containing no relevant scalar operators. 
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1 Introduction 


The conformal bootstrap [1-3], originally successful in elucidating 2D conformal field the¬ 
ories (CFTs), has recently become a powerful method to constrain the operator algebra 
of unitary CFTs also in D > 2 spacetime dimensions. The origin of this development 
is the observation made in [4] that the combined constraints of crossing symmetry and 
unitarity on a 4-point function of scalars can be explored numerically. This method achieved 
an impressive degree of success, for example, by enabling accurate determinations of the 
dimensions of low-lying operators in particular CFTs, such as the 3D Ising [5-8] and critical 
0{N) vector [9, 10] models. While the original ideas of [4] have been developed in myriad 
ways in subsequent works, the basic objects of study have always been 4-point functions of 
scalar operators. The goal of the present paper is to extend the bootstrap toolbox to study 
4-point functions of fermionic operators in 3D CFTs. 

There are many motivations for bootstrapping fermionic correlators. Fermionic opera¬ 
tors exist in many interesting CFTs, though they do not appear in the operator product 
expansion (OPE) of scalar operators. Therefore, in order to access this sector of the operator 
algebra, one must study correlators of fermionic operators. More generally, it is of great 
interest to apply the conformal bootstrap to 4-point functions of operators with non-zero 
spin. For example, studying the 4-point function of a global symmetry current would lead 
to universal bounds on all CFTs that admit the corresponding symmetry, without making 
any additional assumptions on their operator content. Similarly, bootstrapping the 4-point 
function of the stress-tensor would allow for the most general constraints, since the stress- 
tensor exists in any local CFT by dehnition.^ Our numerical study of fermion correlators 
can be seen as a small step towards implementing the bootstrap for operators with spin, as 
in the examples discussed above. 

In implementing the bootstrap for correlators of fermionic operators, or of operators 
with spin more generally, one faces two difficulties. The hrst is that the explicit form of the 
conformal block decomposition for higher-spin correlators is not known in general for D > 2 
dimensional CFTs. The exception occurs in D = 3, where the conformal blocks for external 
operators with integer spin can be determined by acting with certain differential operators 
on the scalar blocks [17]. As we will see, a similar strategy applies for operators with half¬ 
integer spin. The second difficulty has to do with the proliferation of conformal invariants 
that can appear in higher-spin correlators. In particular, the number of conformal invariants 
grows with the spin of the insertions (see, for instance, [18, 19]). As we will explain, a 
4-point function of primary operators of spin-1/2 generally depends on 16 independent 
conformal invariants. Imposing parity symmetry reduces the number of invariants to 8, and 
this number can be reduced further to 5 if we assume the fermions are identical. In our 
numerical analysis, we will focus for simplicity on this latter case and derive the conformal 
block decomposition of a 4-point function of identical fermionic operators in a 3D CFT. 

With the conformal blocks in hand, we then embark on a systematic study of CFTs 
with a small number of relevant operators. By imposing gaps in the low-lying spectrum 

^In some supersymmetric theories correlators of symmetry currents are related to correlators of scalars. 
The numerical bootstrap has been applied to such cases in [11-16]. 
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of parity-even and/or parity-odd scalar operators, we use the logic originally introduced in 
[4] to derive constraints on the dimensions of the first few scalar operators and on that of 
the fermionic operator whose 4-point function we study. We also hnd a lower bound on the 
coefficient Ct that appears in the two-point function of the canonically-normalized stress 
tensor. 

From these studies we find two exciting results. Firstly, the general bound on the 
dimension of the leading parity-odd scalar possesses a severe discontinuity at a fermion 
dimension of ~ 1.27. This coincides with a kink in the general bound on the leading 
parity-even scalar dimension A^. Based on these features we conjecture the existence of 
a fermionic 3D CFT containing no relevant scalar operators. While we do not know of a 
Lagrangian that would give rise to such a theory, this conjectured “dead-end” CFT would 
furnish a concrete example of self-organized criticality [20, 21] in 3D. 

The second result is that, when we allow a second relevant parity-odd scalar in the 
spectrum with dimension Ag./ = 2 + S for small values of S, the resulting allowed region 
for (A^, Ao-) possesses a sharp kink that appears to precisely coincide with the dimensions 
in the 0{N) Gross-Neveu models at large N. This is natural because in the large-A limit 
one has the expansion A^/ = 2 + 32/(37r^iV) -|- .... By tracking this feature at larger 
values of the gap we reveal information about the small-A Gross-Neveu models, including 
the N = 1 theory which is expected to have A/” = 1 supersymmetry. In addition to this 
sequence of kinks, we observe the emergence of a second discontinuity in the allowed region 
at (A.0, Ag.) (1.078,0.565), which we conjecture could also coincide with a 3D GFT with 
fermionic operators and a large scalar gap. We believe that fully isolating these theories 
will require implementing systems of mixed correlators containing fermions and scalars, but 
based on the results of this study the prospects for learning more about these theories using 
the conformal bootstrap looks very promising. 

The rest of this paper is organized as follows. In Section 2 we start by reviewing the 
embedding space formalism for operators with spin and, using this formalism, derive the 
conformal block decomposition of a 4-point function of identical fermions. In Section 3 
we set up the crossing equations and outline the strategy we will follow in our numerical 
study. Next, in Section 4, we present numerical results for bounds on dimensions of scalar 
operators in theories with fermions that satisfy various assumptions, and also present our 
lower bound on Ct- We end in Section 5 with a discussion of our results. 


2 Embedding Formalism for 3D Spinors 

In order to set up the 3D fermion bootstrap, we need an efficient formalism for keeping track 
of the tensor structures appearing in correlators of fermionic operators. We also need to 
calculate the conformal blocks appearing in the expansion of fermion 4-point functions. Our 
approach will be to use an embedding formalism where the spinorial 3D conformal group 
Sp(4, R) is linearly realized. Similar GFT embedding methods have been developed in [17- 
19, 22-30] and various supersymmetric extensions have been developed recently in [31-38] 
and references therein. Details of our group theory conventions are given in Appendix A. 
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We label the 3D coordinates as with /x = 0,1, 2, and use the Minkowski metric 
in mostly plus signature = diag(—1,1,1). The coordinates transform non- 

linearly under special conformal transformations. It is therefore convenient to introduce a 
different set of coordinates that transform linearly under the action of conformal transfor¬ 
mations. Since the 3D conformal group is isomorphic to SO(3,2), we can relate conformal 
transformations to Lorentz transformations in a 5D spacetime with metric = rjAB = 
diag(—1,1,1,1, —1), where the indices A,B run from 0 to 4. The exact relation between 
the generators of conformal transformations and the SO(3,2) generators can be taken 
to be 


D = - . (2.1) 

Let us denote 5D coordinates that transform linearly under SO(3,2) by capital letters, 
The way to embed the 3D coordinates in 5D space is through the projective null 
cone, which is defined as the space of all points that satisfy the condition X ■ X = 0 and 
are identified up to a rescaling X"^ ~ AX"^. It will be convenient to use lightcone coordinates 
X^ = ± X^, and thus write the 5D coordinates from now on as X = (X^, X"*", X“). 

The exact relation between x^ and X"^ is given by 

Xt^ 

= —, X = X+(x^l,x2), (2.2) 

where = rj^yX^x'^. Note that the parameterization (2.2) obeys X ■ X = 0. 

2.1 Embedding of Scalar Fields 

To find the embedding of fields in 5D spacetime we follow the approach of [26]. Consider 
first a real scalar primary field 0(x). Its transformations under the conformal group are 

x[M^",0(x)] = (x"a^ - x^a")(/.(x), 
x[P'^,0(x)] = -a'^0(x), 

i[K^, (j){x)] = {2x^x''dy — x^d^ -|- 2A^x^) 0 (x), 
i[D, 0(x)] = + A^) 0(x), 

where A^ is the dimension of 0. We can relate 0(x) to a scalar field *h(X) defined on the 
lightcone in 5D as: 


where x is related to X through (2.2). Explicit calculation then shows that *h(X) is a 5D 
Lorentz scalar, i.e. that it transforms under 5D Lorentz transformations as 

i[j“, <I.(X)] =(x‘>-X. . (2,5) 
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if and only if (j){x) is a primary scalar field in 3D with dimension A^. 

Note that the 5D held $(X) dehned in (2.4) is a homogeneons fnnction of X of degree 
—A^. This property together with 5D Lorentz invariance restricts the form of correlation 
functions in embedding space. For instance, the two-point fnnction takes the form 

{4(X.)4(X2)) = ^, Xii^-2X,-X,. (2.6) 

^12 

Using (2.4), we can read off 

{^(xO0fe)> = -g;. x ';, (2.7) 

^12 

as expected. In a nnitary theory, we mnst have c,j) > 0. We will conventionally take = 1 
for scalar operators in this work. 

2.2 Embedding of Spinor Fields 

The above procednre can be applied to primary spinor helds as well. The main difference is 
that snch a held transforms in a spinor representation of the donble cover of SO(2,1), 

which is isomorphic to Sp(2, R). Under the fnll conformal gronp, it transforms as 

i[P^, 'ip^ix)] = -d^'ip^ix) , 

V'“(a:)] = {2x^^x’^d^ - x‘^d>^ + 2A,^a:'^) tl)^{x) + 2ix^{M''^^Tp^J^x ), 
i[D, V’“(a:)] = {x^d^ + A^) ^/>“(x), 

where Here, npper (lower) indices a,/3,..., represent fnndamental 

(anti-fnndamental) Sp(2, R) indices that are raised and lowered with the symplectic form 
Hq ,/3 = —see Appendix A for onr conventions. In 3D, the smallest spinor representation 
is a 2-component Majorana spinor, which is what we will focns on. If we take the 3D 
7 matrices to be real, as we do in Appendix A, a Majorana spinor has real components, 
'0"(a;)* = 'ljj°‘{x). 

In order to efficiently keep track of the 3D spinor indices, it is convenient to introdnce a 
set of anxiliary commnting variables Sa and consider the prodnct 

^p(x, s) = Sa 1 p°‘{x) . (2.9) 

The qnantity '4){x, s) contains the same information as the spinor helds '4)^{x), becanse the 
latter can be recovered throngh ip°'{x) = -^ip{x^s). 

Going to the embedding space, we nse the donble cover of SO(3, 2), which is isomorphic 
to Sp(4, R) with generators = —|[r"^,r'®]. For every 3D spinor held 'if^^^x), we 

wonld like to dehne a 5D spinor held T^(X) on the lightcone (2.2), where npper (lower) 
indices /, J, etc. denote fnndamental (anti-fnndamental) Sp(4, R) indices that are raised and 
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lowered with the symplectic form Qjj = . As in 3D, we can also introduce polarization 

variables Si that help us efficiently keep track of the Sp(4, R) indices, 

^(X,^) = (2.10) 

If we wish, we are free to treat Si as a spurionic 5D held that transforms in the anti¬ 
fundamental of Sp(4, R), and similarly to treat as a position-independent spurionic held 
in 3D that transforms in the anti-fundamental of Sp(2, R) and is invariant under dilatations 
(in other words, Sa is a primary spinor held of vanishing dimension). If we do so, then 
4/(X, S') and 'il!{x,s) become 5D and 3D Lorentz scalars, respectively. Just as in (2.4), one 
can check that the relation 

= ( 2 . 11 ) 

implies that 'h(X, S) is a Lorentz scalar in 5D if and only if s) is a primary held in 3D 
with dimension A^. Since we assumed that Sq, transforms as a dimension-zero primary held, 
then we have that 4/(X, S) is a 5D scalar if and only if is an Sp(2,R) spinor primary 
held of dimension A^. 

To hnish the identihcation between the 3D and 5D spinor helds, we can take^ 


S, = ^ ■ (2.12) 

Using the conventions of Appendix A for the embedding of Sp(2, R) into Sp(4, R), one can 
check explicitly that this relation implies that if Si is an Sp(4, R) anti-fundamental spinor 
in 5D, then is an Sp(2, R) anti-fundamental spinor primary held in 3D with vanishing 
dimension, as desired. Notice that Si satishes the transversality condition 

SiX\j = o, x^j = x^{rAYj, (2.13) 


which is invariant under Sp(4, R) transformations and is consistent with the lightcone 
condition X ■ X = 0. Due to this transversality condition, the 5D spinor held is dehned 
on the lightcone only modulo the shifts T^(X) —)■ \k'^(X) -|- X^j0'^(X), where 0'^(X) is an 
arbitrary spinor on the lightcone. 

^Equivalently, we could have written 4'(X, S) = Si'^^{X) = s)/{X^)X: identified 


¥{X) 


1 /-ai“/3V'^(a:)\ 


such that is a 5D spinor if and only if is a primary spinor field of dimension A^. Then 4'(2f) 

would satisfy the transversality condition X^{X) = 0 and the polarization Sj would only be defined 
modulo shifts Sj —>■ Si + TjX'^ j. One can relate this description to the one presented in the main text by 
taking Sj = SjX'’j and 4'^(X) = X^{X). 

With this in mind, we can relate our formalism to that of [19] for traceless symmetric tensor fields. 
They define a vector that satisfies transverseness Z ■ X = 0 and is defined up to gauge redundancy 
Z ^ Z + XX. Such a vector can be obtained as Z^ = ST'^flS'. Note that the gauge redundancy of S 
gives Z^ —>■ Z^ -b TXT^flS = Z"^ + X^{TQS), which is the correct gauge redundancy for Z. Similarly 
Z ■ X = SXilS = 0 since S is transverse. 
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{X, S) satisfies the homogeneity property 

^(aX, bS) = S), (2.14) 

where a and b are arbitrary and independent. Homogeneity, the transversality condition 
(2.13), and 5D Lorentz invariance restrict the form of embedding space correlation functions 
of tl/(X, S). For example, the only consistent expression for the two-point function is 

(4'(X„S,)>I'(V2,S2)>=jc*X^, (2.15) 

^12 

for some constant c^. Here, we used the notation 

(F 1 X 2 X 3 ... ^„) = SuX2^jX/k ■ ■ ■ , (2.16) 

where X^j is dehned in (2.13) and Qjj = is the Sp(4, R) invariant tensor (see Ap¬ 
pendix A). Using (2.11) and (2.12) in (2.15), we obtain, as expected 

= ic^ ^^2 aI+i ■ (2-17) 

^12 

For Majorana fermions, we have G R, as can easily be seen by using the Majorana 
condition = ^|J°' and the fact that complex conjugation interchanges the order of the 
Grassmann variables. (Recall that we work in a basis where the gamma matrices are real.) 
In this paper, we will take = 1 for all external operators. 

2.3 Embedding of Fields of Higher Spin 

The above discussion generalizes to helds of higher spin in a natural way. In 3D, a spin-f' 
held is a totally symmetric tensor in 2i spinor indices: It corresponds to an 

embedding held that is homogeneous of degree — (A -|- £) in X and also totally 

symmetric in its 5D spinor indices. In index-free notation, we contract all the indices with 
an auxiliary spinor s in 3D or with a transverse auxiliary spinor S in 5D: 

0,{x, s) = , 0^{X, S) = Sj.Sj, ■ ■ ■ . (2.18) 

0£{X, S) is also homogeneous in the variable S, with degree 2i. By the same argument that 
led to ( 2 . 11 ), we must have 


0£{X, S) 


(XH 


-Oe{x,s). 


(2.19) 


With the help of (2.12) and 




a2...oi2e , 


X) 


^_ 

{2i)\ dSa^dSa^ . ..dSa2e 
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0{x,s ), 


( 2 . 20 ) 






one can then reconstruct the correlation functions of ( 9 “i“ 2 ... 02 £ from the corresponding 
formulas in embedding space. 

As an example, the two-point function of O^lX, S) in embedding space is restricted to 
take the form 


(0,(XuSi)0,(X2,S2)}=t^% 


{SlS2)^^ 

vAo+e 

^12 


( 2 . 21 ) 


Here cq is real if Oi is real. If i is an integer, we also have co > 0 in a unitary theory for a 
real operator Oi. 

For future reference, we record that when i is an integer, we could have represented the 
spin-£ operator in terms of a rank-£ traceless symmetric tensor of 50(2,1), namely 
This tensor is related to yjg^ 


(- 1 )'^ ( 2 . 22 ) 
n//^e (no!i...a2£f ] 

2^ '0102 ■ ■ ■ '«2^-l“2^'^ ■ 

It is straightforward to show that for spin-1 operators, the 2-point function in (2.21) can 
also be written in the more familiar form 

{C>'‘(ii)0''(^2)> = flJ^. (2.23) 

where — 2x^x'^/x^. The analogous expression for spin-2 

i {i^<^{xi2)r^{xi2) + i^^{xi2)rp{x,2)) - 


For a conserved current, take O = J with Aj = 2 in (2.23), and for 
O = T with Xt = 3 in (2.24). 


{ 0 '“'( 2 :.)C>'”( 2 : 2 )> = ^ 


operators is 

1 

(2.24) 
the stress tensor take 




2.4 Three- and Four-Point Functions 


3-point functions between two scalars and a spin-£ operator of dimension A are uniquely 
constrained up to an overall coefficient to be of the form 


(4>i(Ai)4)2(X2)0,(^3,^3)) = hX20- 


( 53 X 1 A 253 )' 


Ai+A 2-A+.^ A2-Ai-fA-|-^ Ai-A2-|-A+^ 


A 


12 


A, 


23 


A. 


(2.25) 


31 


This form follows from homogeneity of degree —Ai, —A 2 , —(A -|- ^) in Ai, A 2 , and A 3 , 
respectively, homogeneity of degree 2£ in transversality of ^3 with respect to A 3 , and 
5p(4, R) invariance. In addition, it is useful to note that {Ai, A 2 } = (2Ai ■ A 2 )l 4 , which. 
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together with Xk ■ Xk = 0 , restricts the choice of quantities that can appear in between 
(S' 3 (- ■ ■ )S' 3 ) to what is written in (2.25). Note that 

($i(Xi)$2(^2)C?r(^3, ^ 3 )) = (-l)'(<hi(X2)$2(Xi)0,(X3, ^ 3 )), (2.26) 


so in the case of identical operators $1 = $ 2 , fhe 3-point function necessarily vanishes 
if £ is odd. Furthermore, for real scalar operators we have (<Fi(Xi)<h 2 (X 2 )C>£(X 3 , S' 3 ))* = 
(<hi(Xi)<F 2 (X 2 )( 9 £(X 3 , Fs)), so the structure constants A</,^ 02 O appearing in (2.25) are real. 

When we come to 3-point functions containing fermions, we have the new complication 
that multiple tensor structures can appear. In general, we have 


(^i(Xi,5i)^2(X2,52)0^(X3,53)) 




i‘20 


Ai+A2-A-.£+1 A2-Ai+A+.^ Ai-A2+A+£ * 


X 


12 


X 


23 


X. 


31 


(2.27) 


where the index a runs over all possible 3-point structures ra, to be given shortly. These 
structures can be divided into those that are even under parity Xk —)• —Xk and those that 
are odd. A basis for the parity-even structures is given by 


{SiS2){S3X,X2S3y 

yi ’ 

^12 

(5i53)(^253)(^3XiX2^3)'-' 

yt-\ 


(2.28) 

(2.29) 


while a basis for the parity-odd structures is given by 


r-i = 

Ti = 


{S,X^X2S,Y^^^ [X23(5i^3)(^2Xi53) +X33(^2^3)(^1^2^3)] , 

V 2 V 2 

^12 "^23 ^31 

{S^X^X^Sy^ ' [X 23 (^ 1 ^ 3 )(^ 2 ^ 1 ^ 3 ) -^ 13 (^ 2 ^ 3 )(^iX 2 ^ 3 )] • 

'Y^~^ 2 V 2 V 2 

^12 ^23 ^31 


(2.30) 


The dependence on (Xj, Si) in (2.27)-(2.30) can be derived from the same reasoning as that 
presented after Eq. (2.25). The new ingredients here are the transversality of Si and S 2 
with respect to Xi and X 2 , respectively, as well as the Fierz identities 


(^iX2^3){^2^1^3) = -{^1^2)(^3^iX2^3) + 2 {^i^ 3 ) (^ 2 ^ 3 )Xi ■ X 2 , 

( 51 X 352 ){^3XiX2F3) = -2(5i53)(52Xi53)X2 ■ X 3 - 2 ( 5253 )( 5 iX 253 )Xi ■ X 3 , ^ ^ 

which can be used to show that 3-point structures proportional to the left-hand sides of 
these equalities would be redundant. 

When we restrict to the case of identical fermions Ti = ^2 = T, anti-symmetry under 
10 2 places further restrictions on which correlators can be nonvanishing. In particular, if 
£ is even then we have the constraint = 0 , while if £ is odd then = X^^q+ = 

~ (The ± here only serves as a reminder of the parity of the operators O 
that contribute to each structure in (2.27).) In other words, even-spin operators have two 
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structures of even parity and one structure of odd parity, while odd-spin operators have a 
single parity-odd structure. Further, for any spin the Grassmann nature of fermions requires 
(^(Xi,^i)^(X2,^2)0£(X3,^3))* = -(^(Xi,^i)^(X 2,^2)0£(X3,^3)), implying that all 3 - 
point coefficients i^ust be pure imaginary. 

Let us now consider the general structure of 4 -point functions. Scalar 4 -point functions 
are constrained by conformal symmetry to have the form 


($i(Xi)$2(X2)4>3(X3)4>4(X4)) 


A 3 —A 4 Aj^—A 2 

XiA 2 /X24\ ^ giu^v) 

Y I \ W I A 1 +A 2 A 3 +A 4 ) 

-^ 13 / \^ 14 / w-2- W-2- 

w 1 9 94 


( 2 . 32 ) 


where g{u,v) is an arbitrary function of the conformal cross ratios u = and v = 

V 23 X 14 
V 13 X 24 • 

Similarly, fermion 4 -point functions must take the general form 


^l)^2(X2, ^2)^3(X3, ^3)^4(X4, ^ 4 )) 

/X 24 \^ Y.iti9\u,v) ( 2 . 33 ) 

1 ; 1 ) A4+A2+1 A3+A4+1 ) 

\A13y V-^ 14 / W-2--2- 

^12 "^34 

where the tj are a basis of tensor structures that can appear in the 4 -point function. There 
are many choices of bases, but one convenient way to organize the structures is in terms of 
their properties under various exchanges of the coordinates. 

In general, we hnd that there are 8 independent structures of even parity that may 
appear. Anticipating applications to the bootstrap, we will choose 4 structures to be 
symmetric under the exchange 1 <H- 3 , and 4 to be anti-symmetric. A basis for the symmetric 
structures is: 

{S,Ss){S2[X,,Xs]S,) {S2S,){S,[X2,X,]Ss) 

^ 2 Xi ■ X 3 2 X 2 ■ X 4 

{S,X2Ss){S2X,S^) (^1X453) {^2X1^4) {S1X2S3) {82X38^) {S,X,S^) {82X38,) 

^ Xi ■ X2 Xi ■ X4 X2 ■ X3 X3 ■ X4 

{^iX 2 ^ 3 )(^ 2 Xi^ 4 ) {^1X4^3) (^2X1^4) (^1X2^3) (^2X3^4) (A1X4A3) {^2X3^4) 

Xi-X2 ^ Xi-X 4 X 2 -X 3 X 3 -X 4 

(^ 1 ^ 3 ){>S 2 [Xi,X 3 ]^ 4 ) (^ 2 >g 4 )(>gl[X 2 ,X 4 ]^ 3 ) 

^ 2X1 ■ X3 2X2 ■ X4 

( 2 . 34 ) 
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and a basis for the anti-symmetric structures is: 


h = {SiS3){S2S,), 

{S,[X2,X4Ss){S2[X,,Xs]S,) 

^ 4(Xi-X3)(X2-X4) 

{S,X2Ss){S2X,S,) {S^X,Ss){S2X,S,) {S,X2S,){S2X:,S,) {S,X,S,) {82X^8,) 

^ Xi-Xa ^ X^-X^ ^ Xs-Xg ^ X3-X4 

{^lX2^3)(>g2Xi^4) (>glX4^3)(>g2Xi^4) (^1X2^3) (>^2X3^4) (>^1X4^3) (>^2X3 ^4) 

® Xi ■ X2 Xi ■ X4 X2 ■ X3 X3 ■ X4 

(2.35) 


When we restrict to the case that all fermions are identical, \hi = 4^2 = 4/3 = 4^4 = \h, 
there are additional constraints on the allowed structures coming from exchange symmetries. 
Some of these are highly nontrivial and lead to the bootstrap conditions discussed in the 
next section. However, there are also trivial constraints on the allowed tensor structures 
coming from exchanges that leave the cross-ratios u and v invariant: {1,2} o {3,4}, 
{1,3} O {2,4}, and {1,2} -H- {4,3}. Symmetry under these exchanges then forces 

g^ = g^ = g^ = 0. (2.36) 

In other words, restricting to identical fermions means that there are only 5 allowed tensor 
structures. 


2.5 Conformal Blocks 

Now we would like to understand how the 4-point functions described in the previous section 
can be decomposed into conformal blocks, which sum up the contributions of all descendants 
of a given primary operator appearing in the 4/ x 4/ OPE. There are many approaches to 
computing conformal blocks in H > 2, including direct summation [2, 39-43], solving the 
Casimir differential equation [5, 44-47], pole expansions [7, 9], and evaluating monodromy- 
projected conformal integrals [27, 38, 48]. We will here adopt the latter formulation, since 
it will allow us to express the fermion conformal blocks in terms of derivatives of known 
scalar conformal blocks, similar to the approach of [17, 30]. 

Let us briefly review the conformal block decomposition of a four-point function of scalars 
(0i02</>304)- Performing the s-channel OPE, one can write the function g{u,v) appearing in 
the scalar 4-point function (2.32) as a sum of conformal blocks:^ 

g{u, v) = '^ A</,i02oA</,304ci5fA/;Ai2,A34(w, v) , (2.37) 

o 

^Note that the correctness of this formula depends, crucially, on the normalization of the function 
g. In terms of the coordinates r and 6 introduced in [47], Eq. (2.37) holds provided that g{u,v) ~ 
( 1 / 2 P (~l)^(4r)^Pe(cos 6) as r —)■ 0, with bz normalized as in (2.6) and O normalized as in (2.21) with 
CO = 1 . 
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where the sum runs only over primary operators belonging to both the x 02 and 03 x 04 
OPEs. As described in [27], the conformal block of O can be obtained from the integral 


Ai — A 2 

^24 A ^ 5'A,r;Ai2,A34(^;'^) 

V I ^ 1+^2 -^ 3+^4 

^ 12 ’ AV’ (2.38) 

= ^J £''‘A„{®i(A',)<*2(A2)0,(A'„)){6((X„)®s(A'3)<*4(A4))|^, 

where Oi is the shadow operator of dimension 3 — A whose indices are contracted with those 
of Oi, denotes a monodromy projection, and A/0 is a normalization factor. 

Similarly, performing the s-channel OPE in the fermion 4-point function (2.33), one can 
write 


A 


(j)l4>20^(j>3Cj)40 


X 


14 


X 


13 


-^yiy20'^y3V’40 5 'a,£;Ai2,A34(“’ '^) > 

O a,b 


(2.39) 


where the index / runs over 4-point structures, while a, b run over 3-point structures. 
Similarly to the scalar case, the outer sum in (2.39) runs over the conformal primaries 
O that belong to both the 0i x 02 and 03 x 'ijj^ OPEs. The inner sum in (2.39) is new in 
the fermion case; it is present because, for any O, there are several OPE coefficients that 
need to be specihed, as in (2.27). In analogy with (2.38), the conformal blocks appearing 
in the fermion 4-point function (2.33) can be expressed as 


A 


A, 


'lpl1p20'^''lp3llj40 I 


a.b 


X 


14 


"13 


X. 


24 


X 


14 


^/£/5Ui2,A34(“’^) 

A1+A2+I A3+A4+I 


^12 ^ 


^34 ^ 


1 


= ^ / D3Xo(Ti(Xi,5i)T2(X2,^2)O0Xo)) 

^fo j 

X ((P0Xo)T3(X3,^3)^4(X4,^4)) 


(2.40) 


M 


where the index / runs over 4-point function tensor structures. Thus, it is clear that if each 
structure appearing in the 3-point functions {'^i^20i) in (2.29)-(2.30) can be written as 
derivatives of the scalar 3-point functions {^i^20i) in (2.25), then the fermion conformal 
blocks can be computed from the known scalar blocks. This is the approach that we take 
in this paper. 

Concretely, the parity-even structures in Eq. (2.29) can be generated by applying certain 
linear differential operators to {^i^20i). In constructing these operators, we dehne 


6Xk ~ dX^ ’ 


(2.41) 
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and note a few useful identities: 


{S,J-S2){S3X,X2S3) = -2{S2Ss){S,X2S3) , (2.42) 

oAi 

(^2^^i)(^3^1^2^3) = 2(^i^3)(M^3) , (2.43) 

0 X 2 

{S^J--Xs,){S3X,X2S3) = 8 {S^Ss){S2Ss) , (2.44) 

dX\ 0 X 2 

{S,Xs^)[X, ■ Xu) = {S,XuS2) = -{S2Xs,){X2 ■ X,), (2.45) 

oXi 0 X 2 

{SiXXsj)(A'i . X 2 ) = b{SA). (2.46) 


Note that in order for differential operators on Xu, Su to be well-defined, they must 
preserve the ideal generated by the relations Xl = 0,XkSk = 0, {SuSu) = 0. This is indeed 
true for the operators above, though the derivative 7 ^ is not well-defined on its own. 

It can be checked that the 3-point structures appearing in the {^i^ 20 i) 3-point function 
function can be written in terms of the structure appearing in the 3-point function of two 
scalars and a spin-£ operator. Explicitly, we have 


Ai+A2-A-^+1 A2-Ai+A+.£ Ai-A2+A+^ 


X 


12 


X, 


23 


X. 


31 


Vn 


( 53 X 1 X 253 )^ 

Ai-|-A2-A+£ A2-Ai+A+.^ Ai-A2+A+^ 

V 2 V 2 V 2 

_^12 ^23 ^31 


(2.47) 


where we dehned the differential operators 
V, = (5i^2)ni,i, 

^ _ 1 _ LQ\n , (A + Ai + A2-£-4)(A-Ai-A2 

^ 4£(A-l)'‘^^hXi5X2‘^^^ 4£(A-1) 

P 3 = —^ [(^i^^ 2 )n_i 1 - (52^5i)ni _1 

2(A -1) [ 5 X 1 2’2 ^X 2 ” " 


£+ 1 ) 






(s,-hs,)n_, 1 + (Sj As,)n, _i 

6 X 1 2’2 5 X 2 2- 2 


2 ’ 2 

Ai — A^ 


-T>, 


(2.48) 


and Ila^f, applies a shift to the operator dimensions as {Ai, A 2 } —)■ {Ai -|- a, A 2 -|- 6 }. 

Note that Vi and V 2 generate the parity-even 3-point structures, while and 
generate the parity-odd ones. In addition, the operators Vi, V 2 , and V 3 are antisymmetric 
under the exchange 10 2, while P 4 is symmetric. Together with (2.26), these symmetry 
properties imply that in the case of identical fermions Ti = 4 / 2 , we obtain three-point 
functions that obey the anti-symmetry requirement in Xi and X 2 provided that i is even 
when we use "Di, V 2 , and P 3 and that £ is odd when we use P 4 , in agreement with the 
discussion following Eq. (2.31). 
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Defining 


Va = Va\ 


1^3 ,2^-4 ’ 


(2.49) 


the fermion conformal blocks (2.40) are given in terms of the known scalar blocks (2.38) 
according to the prescription: 


^3-^4 




13 


X 


24 


X 


14 


2 


I I\ab 

^l9A,e;Ai2, 


^34 ' 




A1+A2+1 A3+A4+1 


X 




12 

Xi4 


X. 


X 


13 


34 

^ 3-^4 

2 


X 


24 


X 


14 


5'A,r;Ai2,A34 (u, n) 
^1+4^2 *^3+^4 

X 2 ^ 


^ 12 ^ 


"34 


(2.50) 


The explicit formulas for in terms of g are rather complicated, and we will not reproduce 
them here. Writing the fermionic blocks as derivatives of scalar blocks is useful for numerical 
applications. Computing derivatives of scalar blocks is straightforward, for example we use 
the pole expansion derived in [7] to compute the expansion in radial coordinates [47] to 
order Derivatives of fermionic blocks are then obtained as a linear transformation on 
derivatives of scalar blocks. 


3 3D Fermion Bootstrap 


Let us return to the 4-point function of identical Majorana fermions in a parity preserving 
3D CFT. Using (2.12) in (2.33), we can write this 4-point function as 

{'Ipixi, Si)i;{x 2 , S 2 )ij{x 3 , S3)'lp{x4, S4)) = 2A,,+1^ 2A,,+1 Yl v) , (3.1) 

"''12 "''34 I 


where tj = ti{xi, Si) are the 5 different tensor structures that can appear. Crossing symmetry 
under 1 -f-)- 3 gives a constraint 

v^^+^'^tig\u,v) = , (3.2) 

I I 


where the minus sign on the right-hand side comes from the Grassmann nature of fermions. 
In general = Mjtj is related by some matrix M, but in the previous section we have 

chosen a basis of 4-point structures such that | and ti_ | . In this 

basis the crossing relation becomes 






(3.3) 
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or, isolating each tensor structure, 


0 = v) + {v, u ), 

0 = (m, v) — (w, u). 


(3.4) 


Now, the functions g^^{u,v) have a conformal block decomposition: 
g^^{u,v)= 

C>+,£even 0~,leven 

a,b=l,2 ^3 

+ Y1 > 

O-Jodd 

where 0± has parity ±, and we have chosen a basis of parity-odd 3-point structures such 
that the a = 3 structure only allows even spins and the a = 4 structure only allows odd 
spins, as before. For brevity, we have written instead of and will henceforth 

continue to do so. Thus, we can write the crossing equations as 

0= ^ («,,,)+ 5; (a|-)"4a((«.»')+ E (Ay)4'4':i/(«.0. 

0^,1 even 0~,£even 0~,£odd 

a, 6 = 1,2 

(3.6) 


where '“) ± “)■ 

This is the starting point for the usual bootstrap logic. We can exclude assumptions on 
the spectrum by applying a linear combination of functionals af. 


0 ~ ^0+^0+^I± 

I± C>+,£even 
a,6=1,2 

+ (^3tA,£(“W)) + Y (^4W(“’^)) 

0~ ,leven C>“,£odd 


(3,7) 


where we look for functionals that satisfy the constraints 
- Y (^a6Ao(“W)) > 0, 


a,6=1,2 


«/± (^atA,£(“W)) ^ 0, 

(^3W(«W)) > 0 , 


for all parity-even operators with £ even 
for all parity-odd operators with ^ even 
for all parity-odd operators with £ odd. 


(3.8) 
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Recall that in our conventions, all are pure imaginary — hence the extra sign in the hrst 
line above compared to the usual conditions for scalars. The OPE coefficients of the unit 
operator are given by We search for functionals satisfying these constraints by 

approximating the search as a semidehnite program and implementing it in the solver SDPB 
[8]. Details of this implementation are given in Appendix C. 

4 Results 

We can now use the formalism derived in the previous section to derive constraints on the 
space of CFTs. In particular, we consider CFTs with a Majorana fermion 'ip, and focus on 
scalar operators appearing in the p) x pj OPE. We assume a parity symmetry, so that we 
can distinguish between parity-odd scalars, which we denote by a, a', a '\... (in increasing 
order of their dimensions), and parity-even scalars, which we denote by e, e', e",... (also in 
increasing order of their dimensions). 

Using the methods described in Section 3, we hrst derive general bounds on the dimen¬ 
sions of these operators, observing sharp discontinuities that we conjecture to coincide with 
a 3D CFT containing no relevant scalar operators. We then study the consequences of 
imposing gaps in the scalar spectrum, making direct contact with the Gross-Neveu models 
(described below) at large N. Finally we study bounds on the coefficient Ct appearing in 
the two-point function of the canonically-normalized stress tensor. 

4.1 Examples of Fermionic Theories 

While presenting our numerical results, it is useful to keep in mind a few simple CFTs that 
have fermionic operators: 

• Free Theory. The theory of a free Majorana fermion pj has Lagrangian 

^ = ( 4 . 1 ) 

where pj = is the conjugate spinor. The fermionic operator pj has dimension 

= 1. There are no parity-even scalar operators appearing in the 'ijj x pj OPE. 
The only parity-odd scalar appearing in x -0 is 00, which has dimension 2. All 
correlation functions in this theory can be computed via Wick contractions using the 
free fermion propagator 

mxy)Mx2)) a 000 . (4,2) 

• Mean Field Theory. Mean Field Theory is a generalization of the free theory that in 
general does not have a local Lagrangian description. Its operators consist of normal- 
ordered products of a fermionic operator 0 and its derivatives, except that in this case 
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all correlation functions are computed from Wick contractions using the generalized 
free field propagator 

{'iP^{xi)^ 0 {x 2 )) oc • (4.3) 

Mean Field Theory is not properly a local QFT because it doesn’t have a stress tensor. 
However, it satishes the properties of unitarity and conformal symmetry that we study 
in this work. In the ip x ijj OPE there are now parity-even scalar operators with 
dimensions 2A^ + 1, 2Ap -|- 3, 2A^ -|- 5, ..., and parity-odd scalars with dimensions 
2Ap, 2A^ + 2, 2 A 1 ;, + 4, .... In the limit A.^ —)■ 1, we recover a free theory, plus 

additional operators proportional to whose OPE coefficients in the ip x ip OPE 

vanish. As in the free theory, the 4-point function of ip in Mean Field Theory satisfies 
crossing symmetry. 

• Gross-Neveu(-Yukawa) model. Another 3D OFT with fermionic operators is the 
critical point of the Gross-Neveu model [49]. In the Gross-Neveu-Yukawa description, 
one starts with N Majorana fermions ipi (with i = 1,..., A a flavor index) and a 
parity-odd scalar field (p, with the Lagrangian 

1 ^ - 1 1 

C = - ip + g(p)ipi - -d^(pd^(p - -m^cp^ - , (4.4) 

i=l 

where g and A are coupling constants. When N is even, this theory can be studied 
perturbatively in d = 4 — e dimensions (see, for example [50]). It has a critical point 
that can be achieved by appropriately tuning the scalar mass m? (a fermionic mass 
term is forbidden by parity symmetry). This critical point is believed to survive down 
to d = 3, where it can also be studied perturbatively in the 1/N expansion [50-53]. 
Previous contact between the conformal OPE and the large-A expansion of this model 
was made in [54]. 

In the context of this work, we consider a four-point function of the fermionic operator 
ip = ipi- (We leave the study of global symmetries in fermionic GFTs to future work.) 
The dimensions of operators in this GFT are not currently available at finite N. At 
large N, the dimensions of the lowest few operators are shown in Table 1 (see also 
Appendix B). 

• The M = 1 super-Ising model. Another example of a 3D GFT with a Majorana 
fermion ip is the M = 1 supersymmetric Ising model. It is defined as the IR fixed 
point of the UV Lagrangian 

- - ]^d^pd^p - ^pipip - ^ (#2 + hY , (4.5) 

(with the parameter h tuned appropriately), which, when setting = gh/2 and 
A = ( 7 ^/ 8 , is nothing but the N = 1 case of the Gross-Neveu-Yukawa model. This 
Lagrangian has Af = 1 SUSY and can be described in terms of a real superfield 
T, = p + Oip + \d6F with superpotential 

lU = hE + (4.6) 

o 
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Z 2 

0{N) 

A 

0 

— 

1 

1 - 32/{3n^N) + ... 

Ipi 

+ 

V 

1 + 4/(37r2iV) + ... 

V’pV’i) 

— 

Sym\V) 

2 + 32/{37r^N) + ... 


+ 

1 

2 + 32/{37r^N) + ... 

0 ^ 

— 

1 

3 + U/I-k^N) + ... 

(p^ 

(_). 

1 

k + 16k{3k — 5)/{3 ti‘^N) -|- ... 


Table 1: Representations and one-loop dimensions of low-lying operators in the large-3D 
Gross-Neveu models. V denotes the vector representation of 0{N). The dimensions of V'ij 
and 4>^ were computed in [51-53] and reviewed in Appendix B. The dimension of is 

computed in Appendix B. 


Note that a superpotential with quadratic or quartic terms in S is forbidden by parity 
symmetry, since S is parity-odd. In the IR this theory is believed to be described by 
the M = 1 super conformal algebra osp(l|4). 

Just like for the Gross-Neveu-Yukawa CFT, the dimensions of operators in the M = 1 
super-Ising model are not known precisely (but may be experimentally probeable 
[55]).^ They must obey, however, relations imposed by supersymmetry such as = 
Ag- -|- 1/2 = Ag — 1/2 or A^/ = A^-' -|- 1, where a = (j), a' = (j)^, e = 0^, e' = (f)"^ 
are the lowest few scalar operators. The relation A^ = A^. -|- 1 was used in [59] to 
derive A^- > 0.565. This inequality was obtained by intersecting the supersymmetric 
line Ag = Ag- -|- 1 with the bootstrap bounds derived from the crossing symmetry of 
unitary Z 2 -invariant CFTs. 


4.2 Universal Dimension Bounds 

Let us start by computing general upper bounds on the dimensions of scalars appearing in 
the 'ip X 'ip OPE. For the moment, we assume only conformal symmetry, parity symmetry, 
and unitarity. 


4.2.1 The Lowest Dimension Parity Odd Scalar 

In Figure 1, we plot a universal upper bound on A^. (the lowest dimension parity-odd scalar) 
as a function of in any unitary, parity-invariant 3D CFT. The bound starts at the point 
(A. 0 , Ao-) = ( 1 , 2 ), corresponding to the free theory. It then grows monotonically with A.^ 
up to A^ 1.27, at which point a sharp vertical discontinuity occurs, and the bound jumps 
from Ao- ~ 2.9 to A^- ~ 7.7. This striking jump suggests that the value A^ ^ 1.27 has 
special signihcance. We discuss possible interpretations below. 

At the least, we can conclude that any CFT with a fermionic operator of dimension 

"^By contrast, the N = 2 super-Ising model has the exactly known dimension Ag- = 2/3 and contact with 
the 3D W = 2 bootstrap was recently made in [56-58]. 


19 







Upper bound on lowest parity-odd scalar a G 'i/j x 'll; 



Figure 1: Upper bounds on the dimension of the lowest dimension parity-odd scalar 
appearing in the 'tp x ijj OPE, assuming only conformal symmetry, parity symmetry, and 
unitarity. The orange region is allowed, and the white region is disallowed. The black dashed 
line starting at the free theory point (A.^, Aq-) = (1,2) gives the relation among dimensions 
specific to Mean Field Theory, while the dashed line starting at (A^,Ao-) = (1,0.5) gives 
the relation among dimensions expected for AA = 1 SCFTs, assuming ip \s a. superdescendant 
of (T. These bounds are determined using the procedure described in Section 3 (see also 
Appendix C) by performing a binary search in A^- with 10“^ precision. The parameter A 
defined in Appendix C is given by A = 23. 


A .0 < 1.27 must have a relevant parity-odd scalar in the ip x ip OPE. Conversely, a CFT 
with no relevant parity-odd scalars in the ip x ip OPE must have > 1.27. In addition, 
we see that any CFT with a fermion of sufficiently low dimension must have a parity-odd 
scalar in the ip x ip OPE of dimension smaller than 7.7. 


4.2.2 The Lowest Dimension Parity-Even Scalar 

In Figure 2, we show an upper bound on A^ (the lowest dimension parity-even scalar) in 
any unitary, parity-invariant 3D CFT. The bound monotonically increases starting from 
the point (A^, AJ = ( 1 , 3) up to a value of A^ 5.1. At this point, we encounter a change 
in slope which occurs at precisely the same value of as the vertical jump in Figure 1. 

Note that the free fermion theory does not contain a parity-even scalar of dimension 
3, since the only candidate ip$ip vanishes by the equations of motion. However, in Mean 
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Upper bound on lowest parity-even scalar e G x -0 

Ae 



Figure 2: Upper bound on the lowest dimension parity-even scalar appearing in the if) x tp 
OPE, as a function of A^, assuming only conformal symmetry, parity symmetry, and unitarity. 
As —>■ 1, the bound goes to A^ = 3 and has asymptotic behavior A^ — 3 oc (A^ — 1)^/^. 
The bound has a kink at A^ = 1.27, which is the same value of A^ at which the bound for 
parity-odd scalars had a discontinuity, see Figure 1. This bound was computed with A = 23. 


Field Theory we have A,. = 2A.0 -|- 1, and hence there exists a continuous family of unitary 
solutions to crossing symmetry that approach the point (A^, AJ = (1,3). By continuity, 
our bound cannot move below this point, and indeed it attains this optimal value to high 
precision. 


4.2.3 A “Dead End” CFT? 

The kink near (A,^, A^) (1-27, 5.1) in Figure 2 is reminiscent of the kink in scalar dimension 

bounds corresponding to the 3D Ising model [5, 6]. Hence we might guess that there exists a 
3D CFT with a fermion of dimension ^1.27 whose lowest dimension parity-even scalar 
has dimension 5.1. 

The vertical jump in the parity-odd sector is also reminiscent of a feature previously 
encountered in scalar dimension bounds. Specihcally, Figure 1 of [7] shows a sharp vertical 
jump in the bound on as a function of A^. (assuming that A^ saturates its upper bound) 
in 3d CFTs with a Z 2 symmetry. That jump went from Ag.' ~ 2.9 to Ag./ 6.8, and occurred 
for 0.517 < Ao- < 0.52 (at A = 11). At higher values of A, the A^- window shrinks and gives 
the correct value Ag- = 0.518151(6) in the 3D Ising model. The height of the jump also 
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decreases, e.g. to Ag-' ~ 5.4 at A = 19. The correct value of Ao-/ in the 3D Ising model is 
approximately 4.5. 

Reasoning by analogy, Figures 1 and 2 lead us to conjecture that there exists a 3D 
parity-invariant CFT with A^ 1.27 and large anomalous dimensions for both the lowest 
dimension parity-even and parity-odd scalars, perhaps A^ 5.1, and 3 < A^. < 7.7. Note 
that this theory would be a “dead-end” CFT because it has no relevant scalar operators, 
giving an example of self-organized criticality [20, 21]. In particular, it would be completely 
attractive under RG flow, and hence would require no tuning to reach criticality (assuming 
Lorentz-invariance is unbroken).^ We are not aware of a natural candidate Lagrangian for 
this theory.® However, the possibility that we have discovered a new “dead-end” CFT clearly 
merits further study. 

4.3 Imposing Gaps: The 3D Gross-Neveu Models 

With the most general possible assumptions, we have made contact with the free theory, the 
limit of Mean Field Theory as A^ —)■ 1, and a conjectured “dead-end” CFT. Meanwhile, 
the Gross-Neveu models and A/” = 1 SUSY Ising model lie well inside the allowed regions in 
Figures 1 and 2. To see them, we must input more information. 

A natural choice for the Gross-Neveu models would be to organize operators according 
to their 0{N) representations and use the constraints of 0{N) symmetry in the crossing 
equations, as in [9, 10, 61]. We leave this investigation to future work. For now, we adopt 
a simpler procedure: we impose gaps in the operator spectrum and study how the bounds 
change as a function of the gaps. 

Specihcally, we will use a lower bound A^-/ > A“™ as a proxy for N and try to determine 
(A^, Ao-) as a function of A™“. Assuming A^-' saturates its lower bound, the dependence 
of (A^,Ao-) on Ao-/ should be consistent with Table 1 at large N. Note that because we 
are considering a single component 'ip = 'ijji of the 0{N) vector, all representations of 0{N) 
appear in the 'ip x 'ip OPE. In particular, we have a = (p and a' = 'ip(^i'ipj) at large N. 

In Figures 3 and 4, we plot the allowed regions of (A^, Ao-) assuming Ao-/ > A™™ for 
several values of A™“. All allowed regions are consistent with the free theory at (A,;,, Ag-) = 
(1,2). However, the gap in A^i has the effect of carving out the allowed region below the 
free theory, revealing new kinks. The positions of these kinks closely track the large-A 
prediction for the Gross-Neveu models, and hence we conjecture that this family of kinks 
(in the limit A —)■ cx)) interpolates between the 3D Gross-Neveu models. In Figure 5, we 
plot (Ao-', A^) for the kinks in Figure 3, compared with the large-A prediction at 1-loop, 
Ending excellent agreement. 

In Figure 4, we also see that a new kink appears near ^ (1.078,0.565) when 

®Such a theory would also be interesting from the perspective of the AdS/CFT correspondence—an 
AdS 4 holographic dual of this theory would have no tachyonic scalars and hence all moduli would be fully 
stabilized. We thank Eva Silverstein for emphasizing this point. 

^Several examples of 4D dead-end CFTs were constructed in [60]. Their common feature is that they 
are chiral gauge theories, where mass terms are forbidden by gauge-invariance. 
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Allowed (A^, Aq.) assuming A„/ > 2.01,2.03,2.05,2.07,2.09,2.11 



Figure 3: Allowed values of the dimensions (A^,Ao-), assuming A^-' > A™A for G 

{2.01,2.03,2.05,2.07,2.09,2.11}, computed with A = 19. The regions to the right of their 
respective curves (shaded orange) are allowed, while the regions to the left are disallowed. 
The black dashed line shows the relationship between A^ and Aq- using the known 2-loop 
(for Ag-) and 3-loop (for A^) large- A results in Table 1 and Appendix B. The free theory at 
(A^, Act) = (1,2) is always allowed. Below the free theory, there are kinks that closely track 
the dimensions of operators in the Gross-Neveu models at large N. The vertical lines at the 
bottom of the first two curves ensure consistency of the bounds with Mean Field Theory. 


A“™ > 2.3. This new kink is quite robust to changes in A™™. We discuss its possible 
signihcance below. 


4.4 Increasing 

At large N, increasing the gap in Aq-' corresponds to decreasing N. We might hope that 
for big enough A„/, we could obtain information about the theory with N = 1, namely the 
Af = 1 supersymmetric Ising model. In particular, we should identify a feature in the bound 
that coincides with the line predicted by supersymmetry A^ = Ao- + Unfortunately, we 
observe nothing particularly special happening along this line. The kink corresponding 
to larger N Gross-Neveu models becomes somewhat smooth and crosses the SUSY line 
when Act' is slightly smaller than 3, see the upper kink in Figure 7. If this crossover point 
corresponded to the AT = 1 theory, we would obtain the estimate A^ ^ 1.082. However, 
it is possible that using A„r as a proxy for N ceases to work at smaller N. More directly, 
the operator (which we have identihed with a') does not actually exist when N = 1, 
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Allowed (A^, Ag.) assuming A^./ > 2.1, 2.3,2.5, 2.7, 2.9 



Figure 4: Allowed values of the dimensions (A^,Ao-), assuming A^-' > A™/“ for A™“ G 
{2.1, 2.3, 2.5, 2.7, 2.9}, computed with A = 19. The regions to the right of their respective 
curves (shaded orange) are allowed, while the regions to the left are disallowed. The black 
dashed line shows the relationship between A^ and A^- at 2- and 3-loops at large-A. 


so it would be unsurprising if the A > 1 Gross-Neveu kinks are not smoothly connected to 
the N = 1 theory using Ag./ as a proxy for N. The precise fate of the Gross-Neveu kinks at 
small N should become clear when we incorporate the constraints of global symmetry. We 
also expect that the Af = 1 SUSY Ising model will be easier to isolate using a system of 
mixed correlators involving both -0 and a.' We leave both of these investigations to future 
work. 

However, our study of increasing has revealed an interesting feature in the bound 
that appears robust: the kink near (A^, Ag.) ~ (1.078, 0.565) mentioned in the previous 
section. Figure 6 shows the space of allowed dimensions assuming Ao-/ > 3 (equivalently, 
assuming the theory contains exactly one relevant parity-odd scalar), and this kink appears 
prominently. (In fact, it remains present until A^-' ^6.) The zoomed-in Figure 7 makes 
clear that the lower feature is not consistent with supersymmetry, while the upper feature 
(mentioned above) is barely incompatible with supersymmetry for > 3. 

Thus, we are led to conjecture the existence of a non-supersymmetric 3D parity-invariant 

^An exotic possibility is that the W = 1 SUSY Ising model does not actually exist—that the RG flow 
from the free W = 1 theory induced by a superpotential spontaneously breaks SUSY. The theory in the 
IR of this hypothetical flow would contain a free fermion (by Goldstone’s theorem), together with a (possibly 
empty) interacting sector (perhaps a 3D Ising model). We thank Juan Maldacena and Igor Klebanov for 
discussions of this point. 
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Figure 3 Kinks vs. 1-loop large-iV Gross-Neveu predictions 
A.' 



Figure 5: The positions of the kinks in Figure 3 (black points), compared with the 1-loop 
large-prediction A^-/ = — 6 for the 3D Gross-Neveu models in Table 1 (orange line). 

We also indicate the approximate value of N corresponding to each kink. 


CFT with ^ (1.078,0.565) and exactly one relevant parity-odd scalar. The 

proximity of this theory to the SUSY line may suggest that it is closely related to the M = 1 
SUSY Ising model. For example, suppose the Af = 1 SUSY Ising model had a parity-even 
scalar e' that was slightly relevant, < 3. If e' is not the top component of a scalar 

supermultiplet, then deforming the theory by this operator, flowing to the IR, and tuning 
masses appropriately would yield a non-supersymmetric hxed-point with dimensions very 
close to those of the SUSY theory. This possibility could be tested with the bootstrap by 
identifying the M = 1 SUSY Ising theory, determining the dimension and OPE coefficients 
of e', and performing conformal perturbation theory. 


4.5 Central Charge Bounds 


Having explored bounds on operator dimensions, we hnally turn to the “central charge” Ct 
which appears in OPE coefficients of the stress tensor conformal block. We will place a 
general lower bound on Ct as a function of A.^. 

The two-point correlation function of the stress tensor is hxed by conformal invariance 
to take the form (2.24) (with A = 3) up to an overall coefficient. Let us write 




Ct 1 1 

( 471 )^ .2 


(F^(a;i 2 )/""(a;i 2 ) + F" (0:12)/"" (0:12)) 




(4.7) 
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Allowed (A^, Ag.) assuming A^./ > 3 



Figure 6: Allowed values of (A^,Ao-) assuming Ao-' > 3, computed with A = 23. The 
orange shaded region is allowed, while the white region is disallowed. The black dashed line 
shows the SUSY relationship A^ = A^ — The kink in the upper-left corner corresponds to 
the free fermion theory for which (A.^, Aq-) = (1, 2). Figure 7 zooms in on the second feature, 
near (A^,A^) ss (1.078,0.565). 


where denotes the canonically normalized stress tensor, which participates in the Ward 
identity for translations as follows: 

. .. C>„(a:„)) + 5^ ... 0„(x„)) = 0. (4.8) 

i=l ® 

The free boson and free Majorana fermion have = 3/2. 


The Ward identity (4.8) determines the OPE coefficients 
in Eq. (2.27). As we show in Appendix D, we have 


\i _ 
^T,can 


3i(A^ - 1) 


Stt 


X2 

^T,can 


3z 

dvr 


in the 3-point functions 


(4,9) 


In our setup, however, our normalization of operators appearing in the ijj xijj OPE depends 
only on A and £ and is otherwise independent of the details of the OFT we study. In 
particular, we can use the normalization of the blocks explained in Footnote 3, which is 
equivalent to requiring that the 2-point function of the stress tensor (parity-even operator 
with A = 3 and £ = 2) is normalized as in (2.21) or (2.24) with cq = 1. Comparing (2.24) 
with (4.7), we hnd 

(4.10) 
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Allowed (A^, Ag.) assuming A^-/ > 3 (zoom) 



Figure 7: A zoom of Figure 6: allowed values of (A^, Ag-) assuming A^-/ > 3, computed 
with A = 23. The orange shaded region is allowed, while the white region is disallowed. The 
black dashed line shows the SUSY relationship Ao- = A,^ — The feature corresponding to 
the Gross-Neveu models at smaller A^-' has just crossed the SUSY line near A^ ss 1.082. We 
observe nothing remarkable when this happens (and the precise dimensions (A^, Ao-, Ag-') are 
sensitive to A). Notice another feature at ss (1.078,0.565). This feature appears 

already for Ao-' > 2.3 and persists all the way to Ao-' > 6 at the same position in (A^,Ao-) 
plane. 


and consequently 



A 


1 

T,can 


2ti 

y/Cr 


\2 _ \2 27 r 

Arp — Aj. oan ' 


(4.11) 


We can put a lower bound on Ct as follows. In the sum rule (3.7), we isolate the 
contribution of the parity-even spin-2 operator with A = 3: 


\a \ b pU 
^T^T-^ab,3,2 


{u,v) 


a,b=l,2 0+,£even 

a,b=l,2 

0~ ,i even 0~ ,£odd 


where the summation over parity-even operators now excludes the stress energy tensor 
and the identity operator, whose contributions we wrote separately. We now search for a 
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Lower bound on Ct 

CrlCf^ 



Figure 8: Lower bound on the central charge Ct of a unitary CFT containing a fermion with 
dimension A^. As —>■ 1, the bound approaches the free theory value of Ct- The bound 
goes to zero at A^ = 1.47. For larger values of A^ the SDP is infeasible and therefore does 
not produce any bound. This bound was computed with A = 23. 


= 1 , 


functional a such that: 

a,6=1,2 

^ 0; VA > Af, £ odd. (4.13) 

Here, A^ is the lower bound on the dimension of a spin-£ operator, set by unitarity. 
Eqs. (4.12) and (4.11) then imply: 

( 271 )^ 


VA > A^, ^ even. 
VA > Af, i even. 
VA > Af, £ odd. 


< -a/±[E/±o o(M,u)], 


C, 


(4.14) 


where we have used = i5^. Finding a functional a obeying (4.13) places a lower bound 
on Ct- To make the bound as strong as possible, we search for an a satisfying the relations 
(4.13) that minimizes —[F/fg q(m, n)]. This is slightly different from our procedure for 
setting bounds on dimensions, where it was enough just to hnd a functional satisfying certain 
constraints. Nevertheless, the additional task of Ending a functional whose action on a given 
vector is minimal can again be efficiently performed using SDPB. 




Our central charge lower bound as a function of is shown in Figure 8. We normalize 
Ct by dividing by its value in the free fermion theory, Ol?®® = 3/2. The bound has similar 
features to analogous bounds on Ct coming from scalar 4-point functions in four dimensional 
CFTs [62-64], As the fermion dimension approaches its free theory value, —)■ 1, the 
bound on Ct also approaches its free theory value. For larger values of the bound 

becomes stronger, reaching a maximum. In this case, the position of the maximum does 
not coincide with the features observed in the bounds for A^- and A^, and does not seem 
to play an important role as it did in the studies of 3D Ising model. At even greater values 
of A^ the bound goes to zero. After that point, the SDP problem described by (4.13) is 
infeasible, i.e. it is not possible to hnd an a satisfying the constraints in (4.13). Thus, we 
obtain no bound on Ct for those values of A.^, beyond the obvious Ct > 0. 


5 Discussion 

In this work, we set up the 3D fermion bootstrap and explored its numerical implications. 
We hrst developed an embedding space formalism suitable for describing fermionic corre¬ 
lators. We found that conformal blocks for identical spin-1/2 operators are given by the 
action of certain differential operators on conformal blocks for scalars. Using these operators, 
together with differential operators that relate integer spin correlators to scalar blocks [17], 
one can further determine the conformal blocks for any 4-point function in 3D. 

On the numerical side, we have foremost shown that the bootstrap can extract rigorous 
constraints from four-point functions of non-scalar operators. We have obtained general 
bounds on dimensions of low-lying operators in 3D CFTs with fermions and a small number 
of relevant scalars. We also obtained general bounds on the central charge Ct- Our results 
not only provide rigorous constraints on the operator spectrum of CFTs with fermionic 
operators, but also show numerous features reminiscent of those found when applying the 
bootstrap to four-point functions of scalar operators. 

One interesting feature revealed by the fermionic bootstrap is the kink in the parity-even 
bound in Figure 2, coinciding with the apparent decoupling of the leading parity-odd scalar 
in Figure 1. We do not yet know the correct interpretation of this feature, but it is intriguing 
that it may point to the existence of a 3D fermionic CFT with no relevant scalar operators. 
If such a theory is responsible for the kink in Figure 2, it would contain a primary spinor 
operator -i/ of dimension A^ 1.27, and the lowest parity-even scalar appearing in the 
'ip X 'ip OPE would have dimension A^ 5.1. The possibility that these features reveal a 
“dead-end” 3D CFT that gives an example of self-organized criticality merits further study. 

Other interesting features occur when we impose a gap to the second relevant parity- 
odd scalar. By varying its dimension between 2 and 3, we observe a sequence of kinks 
in the (A^, Ao-) plane shown in Figures 3 and 4. When the gap is very close to 2, their 
locations match beautifully onto the dimensions in the Gross-Neveu models at large A, seen 
clearly in Figure 5. At larger values of the gap, we expect that the kink locations make 
precise predictions in small-A Gross-Neveu models. We also observe the appearance of a 
new discontinuity in the allowed region at (A^, A^) ~ (1.078, 0.565), which is robust against 
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making the second parity-odd scalar irrelevant. 

In order to better understand if these discontinuities correspond to specihc CFTs or 
SCFTs, one could pursue three immediate steps:^ 

• One could hope to extend the relation between fermionic conformal blocks and scalar 
conformal blocks to fractional dimensions. By numerically studying the fermionic 
crossing-equations in different dimensions, one could compare the evolution of the 
discontinuities to the results from a perturbative e-expansion (similar to [57, 65]). 

• It is straightforward to extend our analysis to constrain fermionic theories with an 
0{N) global symmetry. As N is varied, we can track the evolution of the bounds 
on operators in each 0{N) representation and again compare with results from the 
large-expansion for the Gross-Neveu models. Such a comparison could help conhrm 
that the kinks in Figure 4 correspond to the hxed-points of the Gross-Neveu model 
with a small number of flavors and in particular determine which kinks in our family 
correspond to integer values of N. 

• Finally, in order to better understand whether theories live at these discontinuities it 
would be fruitful to extend our analysis to mixed four-point functions containing both 
a fermionic operator -0 and a scalar operator 0. This will allow us to impose gaps in 
the fermionic spectrum, opening up the possibility to obtain isolated islands in the 
space of operator dimensions, as was seen for scalar correlators in [7, 10]. We anticipate 
that this analysis will be particularly useful for isolating the Af = 1 super-Ising model. 
E.g., these mixed correlators would allow us to determine the fermionic spectrum in 
the OPE 4> X Ip, enabling us to probe the existence of a conserved supercurrent in the 
spectrum. 

We hope to report on these further investigations in future work. 
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A Group Theory for 3D Spinors 


The 3D Lorentz group SO(2,1) has a double cover which is SL(2, R) ~ Sp(2, R) ~ SU(1,1). 
For us the Sp(2, R) formulation is convenient. It is clear that the smallest irreducible 
representation is a fundamental of Sp(2, R), which has two real components. This describes 
a Majorana fermion in 2+1 dimensions. 

The Lorentz algebra is 


^ - ri'^PMP^) (A.l) 


In the case of SO(2,l), we take the signature to be = diag(—1,1,1) and we have 3 
generators J = ^eabM‘^^ = Ka = where a,6 G {1,2}. 

In the fundamental representation these generators can be written as 

/000\ /00-i\ /0i0\ 

j= 00-i , Ki= i 0 0 0 , A:2=M00 , (A.2) 

\0 i 0 / \-i 0 0 / \000/ 

which satisfy the algebra 

[J, iPi] = tK2, [J, K 2 ] = -tK,, [K,, K 2 ] = -iJ , (A.3) 

and preserve the metric rjM + M'^rj = 0. Here J performs a spatial rotation and Ka 
perform boosts. As usual, the rotation generators are Hermitian while the boost generators 
are anti-Hermitian. 

The fundamental generators of Sp(2, R) (acting on '^") satisfy the same algebra and can 
be written as: 


J = 


1 

2 


0 

—i 




-i 0 \ 

0 i ) ' 



which preserve a symplectic tensor + {Mp^Y'VL = 0, where = 

The (equivalent) anti-fundamental representation (acting on 'ipa = 
with generators J = = J, Ka = VLKa^~^ = —Ka- 


(A.4) 


0 1 

-1 0 


= 

Cq./j'?/’^) transforms 
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The explicit mapping between S0(2,1) and Sp(2,R) is accomplished via a Clifford 
algebra: 


YY + YY = , 

where we can use the explicit real representation 


7° = 


0 1 

-1 0 


7^ = 


0 1 
1 0 


72 = 


The Sp(2, R) fundamental generators are obtained from 


1 0 
0 -1 


(A.5) 


(A. 6 ) 




(A.7) 


Note that in our conventions, the index structure on the 7 ^ matrices dehned in (A. 6 ) is 
( 7 ^)“/ 3 . Indices are lowered by multiplying with from the left, and raised by multiplying 
with from the right (e.g., 7 ^^ = /s and ( 7 ^)“^ = 

The 3D Lorentzian conformal group SO(3,2) has a double cover which is Sp(4, R). 
We would like to identify which Sp(2,R) subgroup corresponds to the Lorentz rotations 
described above. We will write the metric as = diag(—1,1,1,1, —1), where the first 
3 components correspond to SO(2,1) indices. Then the SO(3,2) generators which 

correspond to physical Lorentz rotations and boosts are simply for p, i/ = 0,1, 2. 


The mapping between SO(3, 2) and Sp(4, R) is again realized via a Clifford algebra: 

pApB ^ pBpA ^ , (A. 8 ) 

and spinors transform in a representation of SO(3, 2 ) with generators 


(A.9) 


We can construct a real basis for the matrices explicitly as 


= 


/ 0 

1 

0 

0 \ 


( ° 

1 

0 

0 \ 


/ 1 

0 

0 

0 \ 

-1 

0 

0 

0 

= 

1 

0 

0 

0 


r2 = 

0 

-1 

0 

0 

0 

0 

0 

-1 


0 

0 

0 

1 

7 ^ 

0 

0 

1 

0 

V 0 

0 

1 

0 


VO 

0 

1 

oj 


VO 

0 

0 

- 1 / 

/ 0 

0 

0 

1 \ 


( ° 


0 

0 

1 \ 





0 

0 

-1 

0 

-i-iT 

0 


0 

-1 

0 




(A.l 

0 

-1 

0 

0 

, r" = 

0 


1 

0 

0 




Vi 

0 

0 



V-i 

0 

0 







Y = 


The generators j in the spinor representation satisfy the Sp(4, R) symplectic con¬ 
straint -|- {M^^Y^Vt = 0 with the invariant tensor 


Vtij = = 


/ 

0 

0 

1 

0 \ 


0 

0 

0 

1 


-1 

0 

0 

0 

V 

0 

-1 

0 

oy 


(A.ll) 
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Then in this basis the rotation and boost matrices are block diagonal and are given by: 



/ 

0 

i 

0 

0 \ 


/ -i 

0 

0 

0 \ 

1 


—i 
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0 

0 
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0 

i 

0 

0 

2 


0 

0 

0 

i 

, Ki — — 

2 

0 

0 

i 

0 


V 

0 

0 

—i 

oj 


V 0 

0 

0 

/ 





( ° 

i 

0 

0 \ 

1 

i 

0 

0 

0 

2 

0 

0 

0 

—i 


VO 

0 

—i 

0 / 


A.12) 


In other words, the upper two components of a Sp(4, R) spinor transform like an Sp(2, R) 
fundamental, and the lower two components transform like an Sp(2, R) anti-fundamental: 


T = 



(A.13) 


B Gross-Neveu-Yukawa model at large N 

In this Appendix, we collect known results on the dimensions of low-lying operators in the 
Gross-Neveu-Yukawa model at its conformal hxed point. The Lagrangian of the Gross- 
Neveu-Yukawa model was given in (4.4). At the GFT point, one tunes the mass for the 
scalar held 0 to zero, and one can ignore the quartic scalar interaction as well as the kinetic 
term for <p. After rescaling 0, the Lagrangian takes the form 

1 ^ - 

£ = ^ ■ (BT) 

i=l 

Recall that the Majorana condition in Lorentzian signature is 'ip = In our conven¬ 

tions, 7 ° = icT 2 , so p) = —'tp'^a 2 - 


B.l Dimensions of ?/), 0, and 


The dimensions oiip, 0, and 0^ have been computed in [51-53] at large N. Let us summarize 
some of these results. The dimension of the fermion operator ipi is known up to order 1/N^\ 


— 1 -|- 


+ 


896 32(-668 + Mlvr^ + 3247r2 log 2 - 3402C(3)) 


The dimension of 0 is 

= 1 

The dimension of dP' is 


+ 


2437r6Ar3 


32 32(304 - 277r2) 


+ 0{l/N^ 




277r*m 


, 32 64(632 + 277r2) 

Ai2 = 2 + -^ + 0{1/N^) 


Stt^N 


277r^N^ 


(B.2) 


(B.3) 


(B.4) 
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B.2 Dimension of 


To order 1/iV, the dimension of the operator 0^ can be deduced from the results presented 
above. Indeed, at leading order in 1/iV, the dimension of 0^ equals k. At order 1/N, there 
are only two Feynman diagrams contributing, one scaling as k and one as k{k — l)/2. We 
must therefore have 

= + (B.5) 

for some constants a and b. Comparing with (B.3)-(B.4), we have 

V = ^ + + (b. 6 ) 

For A; = 3, for instance, (B. 6 ) gives 

64 

A^3 = 3 + — + Oil/N ^). (B.7) 


B.3 Dimension of 'ip(^i'ipj) 

We are also interested in the dimension of the 0{N) symmetric traceless operator 
which appears not to have been calculated in the literature. At leading order in N, this 
operator has dimension 2. In the rest of this section, we present the hrst 1/N correction to 
this result, with the combined answer being 

32 1 

(B, 8 ) 


B.3.1 Setup 


To derive (B. 8 ), we hnd it convenient to work in Euclidean signature. The Euclidean 
Lagrangian is the same as (B.l), with the only difference that we should use Euclidean- 
signature gamma matrices, which can be taken to be 70 = < 72 , 71 = cxi, 72 = < 73 . The 
Majorana condition in Euclidean signature is still "0 = —'ijj'^a 2 - 

At N = 00 , the two-point function of '0 is: 




% ~A I Ts 


(B.9) 


In momentum space, this is 


(C(p)^i(-p)) 


^ (7/.*cr2)“V 

" 7 ' ^2 


(B.IO) 
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The effective action for 0 obtained after integrating out the fermions is 




(B.ll) 

with 

1 ^ - 

*j=i 

(B.12) 

Since '0-0 = —' 0 ^cr 2'0 = 

ieQ/ 3 ' 0 “' 0 ^, we have 


4 {A-KY\x-y\ 

(B.13) 

In momentum space. 

V. r ^ N\p\ 

Mp)= , 

(B.14) 

because 

J = -7/2 \p\ . 

(B.15) 

To leading order in 

N we can thus use the propagator for -0: 




(B.16) 

The propagator for 0 is D{p) = l/n(p), or 



D{p) - (0(p)0( p)) - ^ IpI • 

(B.17) 


B.3.2 Anomalous Dimension of 


To compute the anomalous dimension of let us consider the particular case 0{x) = 

The dimension of O is 

/S.Q = 2A^ + T^vertex , (B.18) 


where, in terms of Feynman diagrams, //vertex can be extracted as the coefficient of the 
logarithmic divergence of the vertex correction diagram. Keeping track of all the numerical 
factors and using the propagators (B.16) and (B.17), we have 


hvertexlogA + • • • 

from which we extract //vertex ■ 


1 f d?q tr[ 7 ^ 7 i.]g^g^ 16 8 log A 

2 J (27r)^ N \q\ n‘^N ^ 

8 /( 7 r^A). Using (B.18) and (B.2), we obtain 


(B.19) 




(B.20) 


yielding (B. 8 ). 
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C Implementation in SDPB 


In this Appendix we provide a description of the numerical implementation of the fermionic 
bootstrap using SDPB [8] . In order to implement a semi-dehnite program we limit the space 
of functionals ai^ over which we search over in Section 3, to those taking the form, 


a/* 1/1= E <^!^nST3UU,n 

n<m, 

m-\-n<A 



(C.l) 


with u = zz and v = {1 — z){l — z) and have evaluated the function / at the crossing 
symmetric point z = z = 1/2. 

Applying these functionals to our crossing equation amounts to hnding the {z, z) deriva¬ 
tives of functions appearing in the dehnition of the conformal block (2.40). These 
functions have singularities as z ^ z, the most divergent of them going as {z — z)~^. The 
singularities come from our choice of basis {t/}; the full conformal block is perfectly regular 
at z = z. To avoid dealing with the divergences, we multiply the crossing equation by 
{z — zY before applying the functional a. 

In order to determine the derivatives of the conformal blocks for the fermionic four 
point functions, we have used a Mathematica script to apply the operators Va to the rational 
approximation of the scalar conformal blocks presented in [9]. Thus, the derivatives of the 
fermionic conformal blocks g^^^ can be written as 

« Xf(A)pb-”’’'*(A), (C.2) 

where (A) are polynomials in A and Xe{^) is a positive function for all values of 

A above the unitarity bound. The hat in g should remind us that we actually multiplied 
functions g hj {z — z)^. Consequently, at the crossing symmetric point we can write 




I z=z=l/2 




(C.3) 


where (A) for a, 6 G {1,2} or (a, 6) = (3,3), (a, 6) = (4,4), are linear combinations 

of the polynomials determined in Mathematica using (2.48) and the rational ap¬ 

proximation of the scalar conformal blocks. Using this approximation, we can rewrite (3.8) 
and (4.13) in the form of a polynomial matrix program solvable using SDPB [8], 


Find such that: 

- E ^oo^o,nM,(A„) = i, 

a,6=1,2 

Yab/{A) E 0 for all parity-even operators with i even, 

F 33 ^£(A) > 0 for all parity-odd operators with i even, 

l 44 ,£(A) > 0 for all parity-odd operators with £ odd, (C.4) 
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where the Yab,i are polynomials dehned as 

E (C.5) 

m,n,/± 

for a, 6 G {1, 2} or (a, h) = (3,3), (a, 6 ) = (4,4). In onr applications we take the operator Oq 
on which we normalize to be either the identity operator or the stress-energy tensor. Note 
that becanse of the mnltiplication of crossing eqnation by {z — z)^, some of the constraints 
in (C.4) are identically zero, or their linear combinations are identically zero, i.e. the set of 
constraints is not linearly independent. This can canse instabilities in SDPB, making it rnn 
indehnitely. We want to remove snch “flat directions” and give only linearly independent 
constraints to SDPB. This can be done nnmerically. We can view the set of constraints (C.4) 
as a matrix with rows labeling the constraints and colnmns labeling the components of a 
fnnctional, We then only need to hnd the linearly independent rows of the matrix. 

That can be done for example in Mathematica nsing the bnilt-in RowReduce fnnction. Notice 
that this step needs to be done only once for a given A. 

The fnll description of implementing the polynomial matrix program reqnired to hnd 
can be fonnd in the SDPB mannal [8]. We have used a Mathematica script to manipulate 
the fermionic conformal blocks to obtain the matrix input for SDPB. In order to obtain 
numerically accurate results we have used the parameters presented in Table 2 in our 
SDPB implementation. For A = 19 generating the input hie required by SDPB takes about 
30 minutes (on a single core), while solving each semi-dehnite program takes 25 minutes 
(allowed points) or 100 minutes (disallowed points) on an 8 core machine. For A = 23 
generating the input hie required by SDPB takes about 90 minutes while solving each semi- 
dehnite program takes 3 hours (allowed points) or 14 hours (disallowed points) on an 8 core 
machine. 


D Conformal Ward identities 

In this appendix we study the implications of the Ward identity given in Eq. (4.8) for 
correlators containing fermions. One can multiply (4.8) by a conformal Killing vector 
satisfying oc In a conformal held theory, the fact that the stress tensor is 

symmetric and traceless implies 

f) ^ f) 

— (^i.(x)T'"''(x)C>i(a;i)... On{xn)) -P ^ (5(a; - Xi)iy{xi)^^{Oi{xx)... On{xn)) = 0 . 

i=l * 

(D.l) 

Taking xi = 0, Oi = O, and integrating in x over a small enough sphere of radius e centered 
at the origin, one can extract the integrated OPE 

[ d^nn^Ux)T>^"{x)O{0) = t[Qi, O]{0 ), (D.2) 
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A 

19 

23 

K 

20 

24 

spins 

^19 

S23 

precision 

640 

960 

findPrimalFeasible 

True 

True 

findDualFeasible 

True 

True 

detectPrimalFeasibleJump 

True 

True 

detectDualFeasibleJump 

True 

True 

dualityGapThreshold 

10-25 

10-40 

primalErrorThreshold 

10-25 

lO-ioo 

dualErrorThreshold 

10-25 

0 

1 

0 

initialMatrixScalePrimal (Dp) 

1020 

0 

0 

initialMatrixScaleDual (Dp) 

1020 

0 

4^ 

0 

f easibleCenteringParameter (/^feasible) 

0.1 

0.1 

inf easibleCenteringParameter (Anfeasibie) 

0.3 

0.3 

stepLengthReduction (7) 

0.7 

0.7 

choleskyStabilizeThreshold { 9 ) 

0 

1 

0 

10-40 

maxComplementarity 

lOloo 

10130 


Table 2: Parameters for the computations in this work. Only SDPB parameters that affect 
the numerics (as opposed to parameters like maxThreads and maxRuntime) are included. The 
sets of spins used are S'lg = {0,1, 2,..., 25} U {29,30, 33, 34,37, 38,41,42,45,46,49, 50} and 
523 = {0,1, 2,..., 25} U {29, 30,33, 34, 37,38,41,42, 45,46,49, 50, 59,60}. 


where in deriving the expression we also used Stokes’ theorem and is the conserved 
charge whose associated conserved current is J^{x) = {x). Specializing to Lorentz 

transformations, translations, special conformal transformations, and dilatations, we simply 
replace ^y{x)T^^^[x) with 


{M^p)^,{x) = XpT^ 


flU 


^uTpp , 


iPu)^^ix) = -Tp^ , 

{K^)p{x) = 2xuX^Tpp - x^Tpy , 


Dp{x) = x^Tp^ 


in (D.2), and with M^p, and D, respectively. 


(D.3) 


We are interested in calculating the OPE coefficient between the stress tensor and a 
spinor primary held ^fJ. Using (2.8), Eq. (D.2) becomes 


[ d‘^hh^{M^p)p{eh)i>{0) = - , 

[ cfhh^{P^)p{en)i/j{0) = -^,y^/J{0) , 

Js^ 


/ (fnn^{K^)p{en)'ilj{0)=0, 

Js^ 


/ (fnn^Dp{en)'ijj{0)=A^'ijj{0). 
Js^ 


(D.4) 
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The general form of the OPE x is restricted by the tracelessness and conservation 


of to take the form 


Tf,^{x)^p{ 0 ) = a 


/y-> /Y* 


V’(O) + b 


+ x^xP^pp 


\x\ \x\ 

for some constants a and b. Using the dehnitions (D.3), we have 


^( 0 ) + 


xP{M^,p)p{x)^p{0) = b 


'T^ o/ _ 'y'^ o/ 

^ JaU -lyjyJy Jfjp , , , Q, 


X 


-m+o{x^), 


xP{P^)p{x)^jJ{0) = a‘^^jJ{0) - b ^ ’ 

|x| |x| 

xP{K^)p{x)'ilj{0) = -a‘^'ip{0) - b^^-^p{0) + 0{x), 

|x| |x| 

xPDp{x)ijj{0) = —a^'ipiO) + 0(x°), 

\x\ 


and so 


/ (fnnP{M^p)p{eh)'4!{0) = 

Js^ 


Snb 


lupif^iO) + 0(e) 


e^ / (fnnP{P^)p{en)p{0) = 0(e°), 

Js^ 

e^ [ (fnnP{K^)p{en)'4!{0) = O(e^), 

Js'^ 

/ (fnnPDp{en)p{0) = —87ra'0(O) + 0(e) 

Js^ 


Comparing (D.7) with (D.4), we identify 


® ) b = 


Stt 


IhTT 


The final form of the T x p OPE is 


1- r VpuX^ - ^XpX^ 3 XpXp^p^ + x^xP^pp 

Tp^{x)ij{0) = -1 ,5 V^(0) + 


X 


Ihvr 


X 


</’(•)) + 


(D.6) 


(D.6) 


(D.7) 


(D,8) 


(D.9) 


Let us now compare this expression with what we expect from the 3-point function 
(2.27). For a parity-even operator, we have 


+ A, 


|xi2|'^'^“^+'+l|x23|^+'|x3l|^+' 

2 (xi3)^^"Hx23)'^'""(x31Xi 2X23)"""^ ' ' ' (x3lXi2X23)""^-'""^^ 


O' 


ki 2 r 


2A.0-A+£-l |_ |A+£ |_ |A+£ 


ix 23 r"^ i^sir 


(D.IO) 
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where as usual . From the xz —)■ xi limit of the 3-pt function we can deduce 

the O X 'ip OPE. In this limit, the 3-pt function is 


-““(is)) « Ap- 1 )'-'"^ 


tel) 




P 12 


2A^+1 


ksil 


A+£ 


12 , 


(D.ll) 


+ A^(-ir 


|a;3i| 


A+l 


Using the normalization where ('0“(a;)'^^(O)) = ix^^/ , the OPE contribution of ^|J 

then is 


(t if 

(^ 3 )^/ 3 (^l) _ 


Q;ia 2 


(a^3l) 


a2^-ia2^) 


kail 


A+e 


-'ip^ixi) 


+ i 




(D.12) 


ksp*' 

because this contribution reproduces the 3-pt function in the OPE limit. 

Let us now specialize to the case where O 2 = T is the canonically normalized stress 
tensor. Eq. (D.12) is in this limit 

T"i" 2“3«4(^)^/3(0) ~ + iA|^x^(“ix“2“3^/>"^)(0). (D.13) 


X 


X 


Using (2.22) we can represent the stress tensor as a rank-2 Lorentz tensor: 


Tf,^{x)^pf^{0) ~ -( 07 ^) 

CK1CK2 (^ 7 v) CK 3 Ct 4 




X 


X 


(D.14) 


For the hrst term, we can use 

x'"x^(i(T27;,)aia2(fcr27i.)a3a4(7,x*o-2)^"^""(7p*o-2)“®“'‘^ = [tr( 7 ^ 7 „)tr( 7 ^, 7 ^) -F 2 tr( 7 ^ 7 ^ 7 j, 7 p)] 

= ^ [4:X^Xu - dx^p^i. + Sx^X;,] = -^ {r]^,vX^ - Sx^Xu) . 

(D.15) 

For the second term, we have 

x''xP{ia2'yi,)kik2ii(^2'yu)k3k4{7J(^2y^^^{lpi(^2)^^’'''pJ^^^ 

= ix'^x^ (27<,7^7 p7^V’ + ‘^lalvlp'y^'ip + tr(7p7p)7^7^'0 + tr(7p7i')7<T7M'0)* 

= ^x'^x^ (-27p^7p7^'0 + + 27 pp 7 ^ 7^'0 + (p 77 z/))* ^ 

= ^ [-Ax^r]^^tlj + 6(xpx'^7,4i/ + + I2xf.Xui’Y 

■ 2 2 r 

= -^{x Sx^Xu)^^ + (xpx‘^7,41/ + Xyx’"'^„^)'tp 
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So: 


n ry-> /y»2 rtr* 

T,Ax)m ~ "’“ m ■ ( d . it ) 


\x\ 


4' 


\x\ 


We can compare (D.5) to (D.17) to obtain 


= 


3i(A^ - 1) 


Stt 


\2 _ 

47r ■ 


(D.18) 
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